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In this work we study the cosmology of the general f{T) gravity theory. We express the modified 
Einstein equations using covariant quantities, and derive the gauge-invariant perturbation equations 
in covariant form. We consider a specific choice of f{T), designed to explain the observed late-time 
accelerating cosmic expansion without including an exotic dark energy component. Our numerical 
solution shows that the extra degree of freedom of such /(T) gravity models generally decays as one 
goes to smaller scales, and consequently its effects on scales such as galaxies and galaxies clusters are 
small. But on large scales, this degree of freedom can produce large deviations from the standard 
ACDM scenario, leading to severe constraints on the f{T) gravity models as an explanation to the 
cosmic acceleration. 
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I. INTRODUCTION 

There has been much recent interest in the so-called 
f(T) gravity theory as an alternative to dark energy 
for explaining the acceleration of the universe [1^17]. 
This theory is a generalisation of the teleparallel gravity 
[18, 19] created by replacing T, the lagrangian of telepar- 
allel gravity, by a function /(T). It uses the curvature- 
free Weitzcnbock connection [20] to define the covariant 
derivative instead of the conventional torsionless Levi- 
Civita connection of general relativity 

Teleparallel gravity (see Ref. [19] for a review and ref- 
erences therein) has a set of four tetrad (or vierbein) 
fields which form the orthogonal bases for the tangent 
space at each point of spacetime. They are the dynami- 
cal variables and play the role of the metric tensor field in 
general relativity. The vierbeins are parallel vector fields, 
which gave the theory the descriptor "teleparallel" . It is 
dynamically equivalent to general relativity and so is not 
really an alternative to it, but a reformulation which al- 
lows for a different interpretation: gravity is not due to 
curvature, but to torsion. 

The generalisation to /(T) gravity can be viewed as an 
phenomenological extension of teleparallel gravity (which 
is the special case /(T) ~ T), inspired by the f{R) gener- 
alization (see Ref. [21] for a review) of general relativity. 
However, it has the advantage over f{R) gravity that its 
field equations are second-order instead of fourth-order 
(although it is know that, even though it leads to fourth- 
order equations, f{R) gravity can be ghost free). Yet, it 
also possesses disadvantages. Although /(i?) gravity is 
probably not the low-energy limit of some fundamental 
theory, it does include models that can be motivated by 
effective field theory. In contrast, /(T) gravity seems at 
this stage to be just an ad hoc generalization. 
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Another serious disadvantage of f{T) gravity that was 
pointed out very recently in Ref. [13, 16] is that it does 
not respect local Lorcntz symmetry. From a theoretical 
perspective this is a rather undesirable feature and ex- 
perimentally there are stringent constraints. A Lorentz- 
violating theory is only attractive if the violations are 
small enough to avoid detection audit leads to some other 
significant pay-off. So far, the only such pay-off that has 
been suggested is that f(T) gravity might provide an 
alternative to conventional dark energy in general rela- 
tivistic cosmology. 

The specific models that have been considered in the 
literature [1-12] are rather special as they are tailored to 
reproduce the late-time accelerated expansion of the uni- 
verse without a cosmological constant. However, to do 
so, a parameter in these models is required to be tuned to 
a very low value, comparable with the observed value of 
the cosmological constant. Thus, given the lack of clear 
theoretical motivation for these models, it is rather ques- 
tionable if this can really be considered to be a resolution 
of the cosmological constant problem. 

Nonetheless, given the attention these models have at- 
tracted recently, it seems worthwhile to address their vi- 
ability as alternatives to general relativity in the field of 
cosmology itself, which was their initial motivation. We 
do so by going beyond background cosmology and con- 
sidering linear perturbations. We will show that these 
models behave very differently from the ACDM model 
on large scales, and are, therefore, very unlikely to be 
suitable alternatives to it. 

Cosmological perturbations in these models have been 
considered recently in Refs. [12, 14, 15] as well. How- 
ever, in these papers the field equations are written with 
only partial derivatives and look quite different from the 
usual Einstein equations. Here, we will present a covari- 
ant version of the field equations of f(T) gravity with 
a clear correspondence to the Einstein's equations. We 
will use them to derive the field equations for a per- 
turbed Friedman-Lemaitre-Robertson-Walker (FLRW) 
universe. As /(T) gravity has different dynamical vari- 
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ables (the tetrad fields) from general relativity or f{R) 
gravity (rank- 2 tensorial metric field), we end up with 
new degree(s) of freedom. This fact has been neglected 
before, for example, in Refs. [12, 15]. Here we will show 
how the new degree of freedom arise in the perturbed 
field equations, and numerically assess its effects on the 
linear-perturbation observables, such as the CMB and 
matter power spectra. 

The paper is arranged as follows: in Sect. II we de- 
rive the field equations for f(T) gravity in its original 
form and rewrite it in the covariant form. In Sect. Ill we 
give a detailed introduction to the method of deriving the 
covariant and gauge-invariant linear perturbation equa- 
tions for f{T) gravity theory, and list those equations. 
We focus on a specific model with a power-law functional 
form for /(T) in Sect. IV, and study its background cos- 
mology and the growth of large-scale structure. We sum- 
marise and conclude in Sect. V. Throughout this work 
we use the metric convention (+, — , — , — ) and set c = 1 
and K = SttG, where G is the gravitational constant. 



II. THE /(T) MODEL AND ITS EQUATIONS 

In this section we give a brief description of the /(T) 
model and a detailed derivation of its field equations. 
In contrast to previous works, which wrote these equa- 
tions in terms of partial derivatives of the tetrads, we 
shall do this by expressing them in terms of the Einstein 
tensor plus relevant covariant derivatives of the vierbcin 
field. This approach makes the equations closely resem- 
ble their counterparts in GR and provides a basis for the 
derivation of the perturbation equations in the and gauge 
covariant formalism, which is our final goal. 

Since f(T) gravity is a simple generalisation of telepar- 
allel gravity theory, we shall briefly introduce the latter 
(for a comprehensive review see [19]). 



A. Ingredients of Teleparallel Gravity 



where Einstein convention of summation has been used. 
Eq. (1) implies that in this model g^^, and are 
all dependent on each other, which is important for the 
derivation of the field equations by variation. 

General relativity is built on the Levi-Civita connec- 
tion of the metric 
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(3) 



where a comma is used to denote a partial derivative 
(^ = d/dx^). This connection has nonzero curvature 
but zero torsion. Teleparallel gravity, or the teleparallel 
interpretation of general relativity, instead makes use of 
the Wcitzenbock connection (tilded to distinguish from 



f.a _ i^Q Q l6 



(4) 



which has a zero curvature but nonzero torsion. The 
torsion tensor reads 



^ 13-/ = ^ 7/3 
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(5) 



The difference between the Levi-Civita and Wcitzenbock 
connections, neither of which is a spacetime tensor, is a 
spacetime tensor, and is known as the contorsion tensor: 



fiu 



1 



= -iT P 



rp p _ rpp \ 



(6) 



It is worth pointing out at this point that, based on the 
definition we have given above, the Weitzenbock connec- 
tion, the torsion tensor and the contorsion tensor are not 
local Lorentz scalars {i.e. they do not remain invariant 
under a local Lorentz transformation in tangent space) 
even though they do not have any tangent space indices.^ 
This is the root of the lack of Lorentz invariance in gen- 
eralized teleparallel theories of gravity. 

The Lagrangian density of teleparallel gravity is given 

by 

Ct^^T, (7) 



where 



In teleparallel gravity we have the vierbein, or tetrad, 
fields, hjj {xP)^ as our dynamical variables; Latin indices 
a,b, - ■ ■ run from to 3 and label tangent-space coordi- 
nates; Greek indices /i, z^, • • • run from to 3 and label 
spacetime coordinates. The hP arc both spacetime vec- 
tors and Lorentz vectors in the tangent space. As the 
former (indexed by /^), they are the dynamical fields of 
gravitation, as the latter (indexed by a), they form an or- 
thonormal basis for the tangent space at each spacetime 
point. 

The metric tensor of spacetime, g^jy, is given by 

g,,u = Vabhi.h'l (1) 

where r]ab = diag(l, — 1,— 1,— 1) is the Minkowski metric 
for the tangent space. From this relation it follows that 



= SP hPh'' = 5^ 



(2) 



rp rppf_Ll/rp 



rppl_LUrp ^ rji prpl/fj. 



(8) 



and h = \J—g is the determinant of h!^ with g being the 
determinant of the metric g^^i,. After adding the mat- 
ted Lagrangian density Cm. , variation with respect to the 
tetrad yields the field equations 



hhPSP''^,^ 
\hhlSPP^Tp, 



87rGe' 



(9) 



^ Teleparallelism assumes the existence of a class of frames where 
the spin connection is zero and in which the Wcitzenbock con- 
nection assume the form given in Eq. (1). We choose to work in 
one of these frame. One could introduce a Lorentz covariant for- 
mulation of the theory at the level of the action, but this would 
only change appearances [16] . 
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with related to the usual energy-momentum tensor 
by e^'^ = v'^'^Q'^^h'^ and 



(Tfl 



and clearly 
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B. Field Equation for /(T) Gravity 

The idea of /(T) gravity is simply to promote the T 
in the Lagrangian to become an arbitrary function of T: 



h 



167rG 



fiT). 



(11) 



The field equations are straightforward generalizations of 
those of standard teleparallel gravity just given above: 

/t [dp {hKSj'P) - hKS^^>^Tp,p\ 

+fTThKS,^PdpT + \hh^J{T) = 8^Ge^ (12) 

where /t = df{T)ldT and /tt = d^f{T)/dT^. 

Obviously, if /(T) = T + A with A a constant, then 
Eq. (12) simply reduces to Eq. (9). 



C. Covariant Version of the Field Equations 

The field equation Eqs. (9) and (12) are written in 
terms of the tetrad and partial derivatives and appear 
very different from Einstein's equations. This makes 
comparison with general relativity rather difBcult. In 
this subsection we show that Eq. (9) can be written in 
terms of the metric only and it then becomes Einstein's 
equation. We also present an equation relating T with 
the Ricci scalar of the metric R. These will make the 
equivalence between teleparallel gravity and general rel- 
ativity clear. On the other hand, the tetrad cannot be 
eliminated completely in favour of the metric in Eq. (12), 
because of the lack of local Lorentz symmetry, but we will 
show that the later can be brought in a form that closely 
resembles Einstein's equation. This form is more suit- 
able for introducing the covariant and gauge invariant 
formalism for cosmological perturbations. 

First, let us note that although in Sect. II A the ten- 
sors were all written in terms of partial derivatives, they 
could be rearranged so that all the partial derivatives are 
replaced with covariant derivatives compatible with the 
metric 5^1,, i.e., where Vaff^u/ = 0. In particular, we 
would have 
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where we have used the fact that F^..^ is symmetric in the 
subscripts /3,7. From this it can be readily checked that 



K 



70: 



(14) 



(15) 



in which the square brackets mean anti-symmetrisation, 
and also 



rpUP 



These relations are useful in the calculation below. 

Next, from the relation between F^^ and F^^ as given 
in Eq. ((>) and the fact that the curvature tensor associ- 
ated with the Weitzenbock connection T^^ vanishes, we 
can write the Riemann tensor for the connection Fs,, as 



[19] 



and the corresponding Ricci tensor and Ricci scalar are 



(7U /J-A 



R — K^^P K" — K^^^P K 

— pu J^^pufi 



2V^ i^r"^ 



(16) 



This last equation implies that the T and R differ only 
by a covariant divergence of a spacctime vector. There- 
fore, the Einstein-Hilbert action and the teleparallel ac- 
tion {i.e. the action constucted with the Lagrangian den- 
sity given in Eq. (7)) will both lead to the same field 
equations and are dynamically equivalent theories. 

We can show this equivalence directly at the level of 
the field equations. With the aid of the equations listed 
above, it can be shown, after some algebraic manipula- 
tions, that 

hhP,G\ = {hhPSp^i) - + hh^.T, (17) 

where G^i/ is the Einstein tensor. Substituting this equa- 
tion into Eq. (9) and rearranging, we obtain Einstein's 
equations. If we do the same for Eq. (12) we get 

.frG,, + \g^, [f{T) - frT] 

+fTTS,^pVPT = SnCe^,. (18) 

Eq. (18) can be taken as the starting point of the /(T) 
modified gravity model, and it has a structure similar to 
the field equation of the f{R) gravity. Note that when 
/(T) = T general relativity is exactly recovered, as ex- 
pected. 



III. THE LINEAR PERTURBATION 
EQUATIONS 

In order to study the evolution of linear perturbations 
in the f(T) gravity, we have to linearise the field equa- 
tions. Usually, this is achieved by writing all quantities 
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in terms of the metric perturbation variables, and for 
this we have to use some metric ansatz. In /(T) grav- 
ity, however, it is the vierbein, rather than the metric, 
that is the fundamental field, and it has 16 rather than 
10 independent components. Usually, the six additional 
components correspond to local Lorcntz symmetry, but, 
as mentioned previously /(T) gravity is not invariant un- 
der that symmetry. Consequently, specifying a metric 
ansatz does not necessarily fix all the tetrad components 
[13], and one needs to specify an ansatz for the tetrad 
itself and derive the metric perturbations thereafter. 

Other subtleties emerge here. For example, one can- 
not write the metric in some familiar gauges {e.g. con- 
formal Newtonian) a priori, as these gauges are obtained 
by gauging away certain degrees of freedom in the met- 
ric fields. However, in /(T) gravity the degrees of free- 
dom are different and the lack of local Lorcntz invariance 
means that we can only gauge away 4 of the 16 compo- 
nents of the tetrad due to the invariance of the action 
under spacetime coordinate transformations. 

We follow Ref. [25] and derive the perturbation equa- 
tions in the 3 -t- 1 formalism, in which all the quantities 
are covariant and gauge invariant. This formalism deals 
with physical quantities directly and does not need to 
make a metric ansatz a priori. It is appropriate to use 
it given that we have derived the field equation Eq. (18) 
in the covariant form. It has proved quite useful in stud- 
ies of perturbation evolution in modified gravity theories 
[26-30]. 



A. Covariant and Gauge Invariant Perturbation 
Equations in General Relativity 

The 3-1-1 decomposition makes spacetime splits of 
physical quantities with respect to the 4-velocity u" of 
an observer. The projection tensor Ha^ is defined as 
Hap = dap — UaUp and can be used to obtain covariant 
tensors perpendicular to u. For example, the covariant 
spatial derivative V of a tensor field T^... J is defined as 

v"T^^.;;; ^ h'^h'^ ■ ■ ■ hJh^^ ■ ■ ■ h^vptiwI (i9) 

The energy-momentum tensor and covariant derivative 
of the 4-velocity are decomposed respectively as 

Ba/3 = TTap + 2g(aU/3) + pUaUj3 - pHaf), (20) 
y aUp = (Jap + ■^aP + ^dHafj + UaAp. (21) 

In the above expressions, iTap is the projected symmet- 
ric trace-free (PSTF) anisotropic stress, qa the heat flux 
vector, p the isotropic pressure, (Jap the PSTF shear ten- 
sor, TjJap = '^[aUp] the vorticity, 6 = V"uq. — 3d/ a (a 
is the mean expansion scale factor) the expansion scalar, 
and Aa = Ua the acceleration; the overdot denotes time 
derivative expressed as </> = u'^V a(l^i brackets mean an- 
tisymmetrisation, and parentheses symmetrization. The 
4-velocity normalization is chosen to be u°'ua — 1. The 



quantities tTq,^, qa, p,p are referred to as dynamical quan- 
tities and Oap-^'cOap, Aa as kinematical quantities. Note 
that the dynamical quantities can also be obtained from 
the energy-momentum tensor Qap through the relations 

p = Qapu'^u'^, 
p=-lH^^Qap, 



qa = H^U^Qf,,, 
TTap = H^H^e^,^ +pHap. 



(22) 



Decomposing the Riemann tensor and making use the 
Einstein equations, after linearisation we obtain five con- 
straint equations [25]: 



= \7''{e^i,u''wap); 



nqa 

Bap 
y^SaP 



2VaO 



^^(JaP^^'^Wap-, 



2" 



V^q^ + (P + P)w^z. 



ap 



'^aP " ' 



(23) 
(24) 
(25) 
(26) 
(27) 



and five propagation equations: 



e+le^ ^i/-Aa + ^{p + 'Sp)^0; (28) 
2 - 1 

'^ap + ■^&(JaP — V(qA^) + £ap + -KTTq^ = 0; (29) 

2 

vj + -Ow - \/ ^aAp] ^ 0; (30) 
o 



— K 

2 



T^ap + ■^0'^aP 



K 

2 



{p+p)aap + "^{aQp) 
Sap + 0£ap - V^S,(„£^)^;uH = 0; (31) 



Bap + 9BaP + V''£,(,e^)p>^ 

+ ^V''^.(„£^)^>'' = 0. (32) 

Here, Sap^iu is the covariant permutation tensor, £ap and 
Bap are respectively the electric and magnetic parts of 
the Weyl tensor Wap^iv^ defined by £ap = u^^u^yVafipv 
and Bap ~ ~\u^u''£^ll''yVpapu- Note that the angle 
bracket denotes taking the trace-free part of a quantity. 

Using the definition of the projected derivatives, it can 
be shown that 



+ fv^UpV^u^ - VpU^V^Up) vx (33) 



for any projected vector field Vp {u'^Vp = 0). In the ab- 
sence of vorticity zu^i, (which is true up to first order in 
perturbation because we are considering the scalar mode 
only), the above equation can be written as 



[Vp,V^]wp = -R 



(34) 
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where R^vpX is the spatial 3-curvature tensor defined in 
analogy to the Riemann curvature tensor in the 4D space- 
time (the minus sign is conventional). We can then define 
the corresponding Ricci scalar of the liyperspace perpen- 
dicular to the 4- velocity in the usual way: R — Rf_iy'^ ■ 
With the Einstein equation it is easy to find 



Rk,2kp 



(35) 



The spatial derivative of i?, rja = iaVai?, is then given 
as 



2a ~ 

and its propagation equation by 



(36) 



20 2 - - 

■na + ^Va^ -^daVaV ■ A- anVa^ ■ q. (37) 

Finally, there are the conservation equations for the 
energy-momentum tensor: 

p+ (p + p)6' + V"g„ = 0, (38) 
ga + ^Oq<. + (p + p)Ac. - V„p + V^^TTap = 0. (39) 

As we are considering a spatially-flat universe," the 
spatial curvature must vanish on large scales and so in 
the background i? = 0. Thus, from Eq. (■:!■")), we obtain 



1 



np. 



(40) 



This is the Friedmann equation in general relativity, and 
the other background equations can be obtained by tak- 
ing the zero-order parts of Eqs. (28, 38), yielding: 



+ Je' + f(p + 3p)=0, 

p + {p + p)e^ 0. 



B. Generalisation to the /(T) Gravity 



(41) 
(42) 



In order to make best use of the formulae obtained for 
general relativity, we can consider the modifications to 
the Einstein equation in /(T) gravity as a new effective 
energy-momentum tensor Q'^JJ in addition to that of the 
fluid matter, 8^^. Eq. (18) can then be rewritten as 



(43) 



in which 



IT 



IT 



2f: 



-9iJ.u [/ - ItT] 



h 



(44) 



See Ref. [31] for a discussion on hyperspherical and hyperbolic 
universes. 



As already mentioned, here we have to work with the 
tetrad and not just the metric, so the setup will be 
slightly different than that of general relativity. 

Since we intend to investigate the perturbation evolu- 
tion in an almost Friedmann universe, let us first consider 
an exact Friedmann universe: there is no special spatial 
direction and the fundamental observer's world line is 
in the time direction. Assuming that the comoving ob- 
server's frame is aligned with the frame defined by the 
tetrad in tangent space we have Hq = u^, where is 
the 4-velocity of the fundamental observer, and the /if 
{i — 1,2,3) arc three orthonormal vectors in the 3-space 
of the fundamental observer (here we use an underline 
to denote components of the Lorcntz index). If we de- 
fine Ua = Ku^, then 1 = g^.uf'u'' = ijabhf.hluf'u'' = 
VabWU''. Note that in this case = Sq. 

In an almost Friedmann universe, the above symme- 
try is at best an approximation, and will not coin- 
cide exactly with it'' but could differ slightly. Instead of 
U" = (1,0), we will have U" = {U^.U'-) where the W- 
are small, and r]abU°'U^ = 1 implies that U- = 1 up to 
first order in perturbation. As U- 



hjiu^, we can write 



aSu, where perturbation vector and 



M''ae„ — 0. As it will turn out, all the information we 



need to know about hj, is that /it. 



U- is first order 



in perturbation and hj^. This suffices to show that, to 
this order of perturbation, a^^ is the only new physical 
degree of freedom with respect to general relativity. It 
was expected to appear due to the lack of local Lorentz 
invariance [13]. Detailed calculations in support of these 
statements, as well as explicit derivations of the pertur- 
bative expressions for quantities entering the field equa- 
tions can be found in the Appendix A. Here, we will 
only quote the results of these calculations. Obviously, 
the only quantities that are not already present in gen- 
eral relativity are T and S^^pVT . Up to first order in 
perturbations, we have 



T : 



2 4 ~ 

3 3 



(45) 



V • a; ) Hp,i, - ^tu^Rf, 



-T (^o-^^ + TUpy + V(,iae^) -I- V[^a3t,] 



(46) 



Here V • a: = Vaj^ and i?^ satisfies V^Rf^ = R. 

Using the definitions given in Eq. (22), it is straight- 
forward to obtain 



1 

Ft 



Ut - iW + ^if- ItT) 



1 



(47) 



1 

h 

2 1 
sIt^ 



[It - l)npf - - (/ - hT) 



JttT 61 + V • a; 



(48) 
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eff 



KTT 



eff 



1 

h 
1 

h 
1 



1 

(/t - - T^hrTRa 



ifr - l)nqi - -IttOVT 



ifr - l)^7Tip - IttT [a^^ + V 



(49) 
(50) 



(51) 

up to first order in perturbation. There are two differ- 
ent expressions for q'^^-^, which is because the quantity 
Si/fip^^'T is not symmetric a priori^ but the field equa- 
tions require its antisymmetric part to vanish. 

We are also interested in the density and pressure per- 
turbations, and these can be obtained by differentiating 
Eqs. (47, 48): 

^ [(1 - fT)fiVaP^ + /ttTVc^t] , 



h 



(52) 



1 

h 

4 

3 



(/t - l)KVa/ + -e'^OfTTT^c.T 



(53) 



Eqs. (49, 50, 51, 52, 53), together with the equations 
given in Sect. Ill A, are all we need to study the pertur- 
bation evolution in /(T) gravity. 



C. Scalar Equations in f{T) Gravity 

Our formalism has so far been as general as possible. 
Now we will focus exclusively on scalar perturbations and 
perform the following harmonic expansions of our pertur- 
bation variables 

Vo.P = Y.-xQl, ^^p^Y.-^'Q'c. 

— ' a — ' a 

k k 
fe k 

k k 

fc k ^ 

^3 



£o.p = -Y.-<^Qip (54) 



in which Q'^ is the eigenfunction of the comoving spatial 
Laplacian a^V^ satisfying 



Qi,Qip are given by = f V„Q^Q^^ = fV(„Q^^. 

In terms of the above harmonic expansion coefficients, 
Eqs. (24, 26, 29, 31, 36, 37) can be rewritten as [25] 

^fc^(CT-Z) = Kga2, (55) 



k^(j) = -^Kfl^ [fc(n + X) + iUq] , 



k{cj' + Her) = k^{(t) + A) - -Ka^n, 



(56) 
(57) 



k^{(j)' +n<p) = -na^ [k{p + p)(T + kq - W - nn] , (58) 

k^7] = kXo^ - 2knZ, (59) 
krj' = -Kqa^ - 2knA (60) 

in which % = a' / a ~ ^aO and a prime denotes the deriva- 
tive with respect to the conformal time t [adr = dt). 
Also, Eq. (39) and the spatial derivative of Eq. (38) be- 
come 

q' + mq+{p + p)kA-kXP + ^kU^O, (61) 

X' + 3h'{p + p) + 3'HiX + XP) + kq = (62) 

We shall always neglect the superscript *°* for the total 
dynamical quantities and add appropriate superscripts 
for individual matter species. Note that 



h' = -kZ - UA. 
o 



(63) 



and p,p,X, XP, q, 11 with superscripts ^ or '^^f are the to- 
tal quantities (fiuid matter plus correction terms). The 
harmonic coefficients X'^^ ^ , X^''^^ ^ ,q''f ^ ,W^^ ^ can be de- 
rived from Eqs. (52, 53, 49, 50, 51) such that 



/tk {Xf + X^^f) = nXfa^ + 2A^kn^{Z + is). 



(64) 



fTK {XP'f + XP-^ff) = KXP'^a'^ 

- ^[8kn im' - H^) {Z + as) 

-f %kV?-{Z + a;)' + 24^^ (H' - U^) A\ 



96 



^-kV? (H' - -H^) (Z + a;), 



(65) 



/tk [q^ + g'-^^) = nq^a^ 



a" 

nq^ c? - 12^kn {W - H^) ri 



lllLk^U^^Z + x) (66) 
. /tt , 

(67) 

fTK{nf + n''^f) ^ kU^ a"^ (68) 

- \2^kn in' - n^) (a + x). 

This completes our derivation of the scalar mode co- 
variant and gauge-invariant perturbation equations for 
/(T) gravity, and we have one extra dynamical degree of 
freedom ae. It is now straightforward to choose a gauge, 
and as an example the perturbation equations in the con- 
formal Newtonian gauge are given in Appendix B. 
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FIG. 1. (Colour online) The background evolution for the /(T) gravity model with /(T) = T - / {-T)". Upper-left 

Panel: the fractional energy densities for matter (f^m), radiation (ilr) and the effective dark energy {Qde = 1 — ^m. ~ ^r), as 
functions of the cosmic scale factor a, which is normalised to 1 today. Upper-right Panel: the total effective equation of state 
Weff = — 1 — sis a function of a. Lower-left Panel: the ratio between the Hubble expansion rates for the /(T) gravity 

model and for the ACDM paradigm, as a function of a. Lower-right Panel: /t — 1 as a function of a. Here, results are shown 
for n — (black solid curve), 0.1 (green dotted curve), —0.1 (cyan dashed curve), 0.2 (purple dash-dotted curve) and —0.2 
(pink dash-triple-dotted curve). Note that n = corresponds to the ACDM paradigm. The relevant physical parameters are 
fim = 0.257, ftr = 8.0331 x 10"^ and Ho = 71.9 km/s/Mpc. 



IV. NUMERICAL RESULTS 

For a quantitative analysis of the evolution of cosmo- 
logical perturbations in f(T) gravity, one needs to con- 
sider a concrete class of models. Since the motivation 
for considering f(T) gravity was based on the sugges- 
tion that it could account for the late time cosmic speed- 
up without the need for dark energy, it makes sense to 
restricts ourselves to models that exhibit this property 
(we have expressed out reservations about the theoreti- 
cal motivation of a general /(T) theory in the Introduc- 
tion). Thus, we focus on the class of models that can be 
parametrized as 



where n is some real number. The /x parameter will be 
fixed to such a value so that the model can reproduce the 
late time accelerated expansion of the universe. The mi- 
nus sign in (— T)" has also be chosen with some foresight, 
as T = ^f^*^ = —6H^ < in background cosmology. 
Our aim is to examine if this particular class of model 
which can reproduce the background cosmological evolu- 
tion of the ACDM model is also compatible with large 
scale structure evolution. Such a Lagrangian has been 
studied previously by [5, 14] but in different contexts. 



2(n+l) 



(69) 



8 



Evolution of e 




k = 0.01 h Mpc"' 



_k = 0-1 h Mpc"' 



FIG. 2. The time-evolution of frame-independent quantity 
e = st + aJor the model with /(T) = T - / {-T)" and 

n = 0.1, on different length scales, characterised respectively 
by k/{h Mpc"^) = 10"* (solid curve), 10"^ (dotted curve), 
10~^ (dashed curve) and 10 ^ (dash-dotted curve). 



A. Background Evolution 

In background cosmology, the modified Friedman 
equation is given as 



3i/^ = K (p 



f 4- p<^ff 



= ^«/-;t^(/-/tT). (70) 

Using the fact that T — ~6H'^, this equation could be 
written as 



2(n+l) 

T-{l + 2n)i^——=2npf, 



(71) 



according to which we could fix fi by assuming that the 
present fractional energy density for "dark energy" is VIa- 



1 



2(l+„) ^ (6iJ 

1 + 2n ^ ° 



2\ 1+" 



(72) 



Here Hq is the present Hubble expansion rate. The mod- 
ified Friedman equation can then take the form 



Ho 
H 



(73) 



Here the second term in the right-hand side represents 
the energy density of an effective dark energy component. 
Given the value for n, we can solve the algebraic equation 
Eq. (73) to find the expansion rate of the Universe at any 
earlier time. 

We have considered five different values for n, with 
n = 0.0, ±0.1, ±0.2, and summarised the results for the 
background evolution in Fig. 1. The upper left panel 
shows the fractional energy densities for matter, radia- 
tion and effective dark energy respectively. The black 



solid curve (n = 0) is the ACDM paradigm. Hq/H in- 
creases until it reaches its current value 1.0, so a positive 
n (green dotted and purple dash-dotted curves; same be- 
low) means the energy density of dark energy was lower 
in the past. The opposite is true for a negative n (cyan 
dashed and pink dash-triple-dotted curves; same below). 
This behaviour is as predicted by Eq. (73). 

For positive values of n the energy density of the "dark 
energy" increases in time, which implies that its pressure- 
density ratio should be less than —1. Given that we nor- 
malise the "dark energy" fractional energy density by its 
value today, at earlier times it will be lower in the /(T) 
gravity model than in ACDM and so the universe will 
expands slower than in the latter. The effect on the to- 
tal effective pressure-density ratio of all matter species, 
which is defined as Weff == — 1 — jjjr, is shown in the 
upper-right panel of Fig. 1. 

Meanwhile, since for positive values of n the dark en- 
ergy (and therefore the total energy) density was lower 
in the past than in the ACDM paradigm, the Hubble ex- 
pansion rate for the former must be lower too, as can be 
seen from the lower- left panel of Fig. 1 . Note that at very 
early times the expansion rates in these two models are 
almost the same, because the effect of the /(T) correction 
(or the cosmological constant) is negligible then. 

Finally, we shall find that the quantity /t = df/dT is 
important in the f{T) gravity model and so have plotted 
its evolution in the lower- right panel of Fig. 1 . Clearly 



/t 



1 



2(n+l) 



(-T) 



n+l 



(74) 



and therefore must be negative for positive values of n, 
and vice versa. Again, at very early times /t — 1 ~ 
because \f{T) — 1| ^ |T|, and the deviation of fx from 
unity only becomes large at late times. 

These results show that as long as \n\ is close enough to 
0; the deviation of the /(T) gravity model from ACDM 
is small but the background expansion rate could provide 
a weak constraint on the model parameter n. 



B. CMB and Large-scale Structure 

Having fixed p. in order to reproduce the desired back- 
ground evolution, we are ready to consider the evolution 
of linear perturbations. These could place much more 
stringent constraints on the model parameters. 

In the section above we gave the covariant and gauge 
invariant linear perturbation equations for general f{T) 
models. In order to solve these equation numerically we 
must specify a gauge (or reference frame). As usual, we 
choose to work in the CDM frame (that is, the reference 
frame of an observer comoving with dark matter fluid), 
which is characterised by vcbm = ^ = 0, where vqum 
is the peculiar velocity of the dark matter fluid and A is 
the acceleration of the observer. 

Next, we need to determine the behaviour of the new 
degree of freedom ae. In most modified gravity theories, 
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FIG. 3. (Colour online) The power spectra for the large-scale structure of the /(T) gravity model with /(T) = T — 
^2(i+n)yj-_/jnyi Upper-left Panel: the CMB spectrum for different values of n - (black solid curve), 0.1 (green dotted 
curve), —0.1 (cyan dashed curve), 0.2 (purple dash-dotted curve) and —0.2 (pink dash-triple-dotted curve). Upper-right Panel: 
the same as the upper-left panel, but for the matter power spectrum at redshift (today). Lower-left Panel: the late-time 
evolution of the dark matter density contrast Acdm on different scales (as indicated besides the curves); three values of n have 
been considered - n = 0.0 (solid curves), 0.1 (dotted curves) and —0.1 (dashed curves). Lower-right Panel: the same as the 
lower-left panel, but for the late-time evolution of the gravitational potential (f) on different scales. The physical parameters 
are the same as listed in the caption of Fig. 1, and three species of massless neutrinos are used. 



this will be governed by a dynamical equation. In the 
/(T) gravity, however, its value is given by a constraint 
equation. This is a consequence of the fact that the right- 
hand side of Eq. (12) is not a priori antisymmetric, but it 
is required to be as a consequence of the field equations. 
This leads to the two different expressions in Eqs. (GG) 
and (G7), which imply that 

kn{z + is) ^ ^ {n' - n^) . (75) 

This equation can be used to determine a; in terms of 
Z, rj and background quantities. 

We can then eliminate cB in all the relevant perturba- 
tion equations. Nonetheless, it is interesting to see how 
the new degree of freedom '<& evolves in time on differ- 
ent length scales, and this is shown in Fig. 2. Since te 
is not a gauge invariant quantity, what we have plotted 



is e = ae -t- (T, which is gauge invariant. Note that in the 
conformal Newtonian gauge, in which ct = (c.f. Ap- 
pendix B), the quantity e coincides with te. We show the 
results for n = 0.1 in Fig. 2. We see that e decreases in 
time, and the decrease becomes more rapid as one moves 
to smaller scales (bigger A:'s). Therefore, we expect any 
deviations from the ACDM model to be more important 
on large scales than on small scales. We will confirm this 
below. 



We can now examine the growth of the dark-matter 
density contrast in the context of the f{T) gravity model. 
For simplicity, we shall assume that the universe is filled 
with dark matter only, which is a fair approximation at 
late times. Taking the spatial derivative of the Raychaud- 
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huri equation one gets [32] 

kz' + mz - k^A 



(76) 



+3 {n' -n')A = -- {X + SA-P) a', 

in which k is the wavenumber and A", include con- 
tributions from both the dark matter and the /(T) cor- 
rections. In the CDM frame A = and the conserva- 
tion equation for dark matter gives A' ~ —kZ, where 
A = Xdm/pbm is the dark matter density contrast. 
Then, Eq. (76) can be rewritten, by manipulating our 
set of perturbation equations, as 



A'' 



(1 - 2C) HA' = 



fr 



A (77) 



with C defined by 



C = 216 



fTTr/a,* 



-216- 



fc2 _ 36 



/t 

y2-T7' 



•H2 (w-n^) 



fTT/a^l56n^n' - 2AHH" - mW - 48-H'' 



Clearly, on very small scales, where k ^ H,H' /H and 
%" 1%^ we have C — > and Eq. (77) reduces to that in 
the ACDM model, only with the value of the gravita- 
tional constant rcscaled by 1/ fx- On very large scales, 
in contrast, we can neglect k"^ in the expression for C, 
and Eq. (77) becomes very complicated, leading to large 
deviations from ACDM. 

One should also be able to derive an evolution equa- 
tion for the gravitational potential 4> defined in Eq. (54) 
(indeed this will be easier if we use the Newtonian gauge 
potentials given in Appendix B), but we shall not do that 
here. 

In Fig. ' wc show some results for the linear perturba- 
tion evolutions in the /(T) model studied here. Clearly 
both the CMB and matter power spectra (for all choices 
of n except for n = which corresponds to ACDM) blow 
up on large angular scales (small t or small k), which is 
consistent with the above analysis that the evolution of 
matter density perturbations (and therefore the gravita- 
tional potential) on large scales is very different from the 
ACDM predictions. On small scales, however, the f{T) 
model gives similar predictions as ACDM, which is as 
expected. 

To see more clearly how the growth of the dark mat- 
ter density contrast and the growth of the gravitational 
potential have been modified, we have plotted them in 
the lower panels of Fig. 3. For Acdm, the difi'erence be- 
tween the /(T) models (with n = ±0.1) and the ACDM 
is within ~ 10% on small scales (fc > O.OOl/i Mpc~^) be- 
cause the effective gravitational constant is rcscaled and 
the cosmic expansion rate is modified as well. But on 



very large scales {k < O.OOOl/i Mpc"-'^), the difference 
becomes very significant. The same happens to 4>- 

These results are expected to remain qualitatively true 
for other choices for the function /(T), if they are made 
so as to explain the late-time acceleration of the universe. 
This can be seen from the expression for C, which shows 
that the large-scale deviation from ACDM is inevitable 
whenever fxT and/or /ttt sltc nonzero. 

The results suggest that /(T) gravity models which 
are proposed as an alternative to dark energy could face 
severe difficulties in being compatible with observations 
regarding large scale evolution. The expectation that 
that linear perturbation analysis gives better constraints 
than the consideration of background cosmology alone is 
clearly confirmed here as well. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have given the modified Einstein equa- 
tions for general f{T) gravity models in a covariant for- 
malism, and derived the covariant and gauge-invariant 
perturbation equations in the 3 + 1 formalism. The per- 
turbation equations take full account of the extra degrees 
of freedom in the /(T) gravity theory (the importance of 
which was first discussed in Rcf. [13]) up to linear order. 
The equations in specific gauges can then be obtained 
straightforwardly as shown in Appendix B. 

For a general /(T) theory it turns out that no new 
degrees of freedom appear at the background level, and 
the modified Fricdmann equation is simply a nonlinear 
algebraic equation in the Hubble rate H that can easily 
be solved numerically. At the linear order in perturbation 
there is a new vector degree of freedom (as a consequence 
of the lack of local Lorentz symmetry, as pointed out in 
Ref. [13]). However, at this order the equations include 
no time derivatives of this vector, which just satisfies a 
constraint equation. 

After developing the general formalism and deriving 
the perturbed equation at linear order, we restricted our 
attention to scalar perturbations. We then considered 
a broad class of /(T) theories which are representative 
examples of models that could account for the late-time 
acceleration of the universe, as proposed in the litera- 
ture. We studied in detail their background cosmology 
and the evolution of linear perturbations. We were able 
to determine the new degree of freedom algebraically in 
terms of other curvature perturbation quantities. We also 
derived the evolution equation for the dark-matter den- 
sity contrast A in a dark-matter-dominated universe, and 
showed that it resembles that of ACDM on small scales, 
but gets significantly modified on large scales. The large- 
scale CMB and matter power spectra blow up, signalling 
a serious viability problem for any f(T) models that are 
able to account for the accelerated expansion of the uni- 
verse at the background level. We have argued that this 
conclusion is robust and holds true for other choices of 
/(T) unless Jtt = /ttt = at late times. 
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Our result clarifies the effects of the new degree of free- 
dom in the f(T) gravity model at the linear perturbation 
level, and we have seen here that only one extra degree 
of freedom arises. An interesting question is whether 
further degrees of freedom will enter into the field equa- 
tions, and if so, whether they are well behaved, when 
followed beyond linear perturbation. This will be inves- 
tigated elsewhere. 



Appendix A: T and Svi^pVT up to First Order 

We give the perturbative expansions and calculations 
needed to derive Eqs. (45) and (4()). First, we need to 
express the covariant derivatives of hj^ and ae^ in terms 
of perturbation quantities in the 3-1-1 formalism. Using 
the definition V ^hr^ = H'^H^V ah-^ it is straightforward 
to show 



(Al) 



up to first order. Note that h\, and V ^hj, are both first 
order. Similarly, for s^, which is itself first order, we 
have 

V^a5,. sa u^sii, + V^a3,. - '^Oui^aif^. (A2) 

Next we consider T. Using Eqs. (14) and ( Hi) we find 
that 

= (V^K) (V^K) - ,M C^^K) i^uK) . (A3) 
Then, given Eq. (Al), we can show that 

to first order, and therefore 



2 A - 

3 3 



(A4) 



Similarly, it can be shown that 

S.^pS/PT - uPfS.^^p + S.ppVPT (A5) 
where, to first order. 



(A6) 



According to Eq. (15), 



with, to the same order, 

Uph^^\/^hP « (61 4- VPa2p ) H^^ - {Ap 4- stp) 



- I CT^iy 4- n7 



-gpuuPh'^.Vxhl^^e + VPs.p) gp,, 



(A7) 



Clearly now we need to calculate ( h-pV ) . Note that 



this is a vector which is first order in perturbation, and 
i^Ji^xhf^ = U^- (y\hi^ « 0, which means that the 

part of (^hjiVxh^^ which is parallel to it^ vanishes up to 
first order, so we need to consider only the part perpen- 
dicular to uf", (h-pS/\h^^ ^ = Tp. 

In order to find an expression for T^, consider Eq. (16), 
rpiyp^ = XA";; and Eq. (14), which leads to 

2V^ (/t^V^/i^) =-R-T. (A8) 
Using Eq. (A 4) and the relation 



-26- -r + 2VP-Ap - R 
o 



(A9) 



(AlO) 



to first order [29], we have that 

2^ (h^^V^K) w 2^ + 26*2 - 2V • A 

. 4 . 
+R+-eV -ee. 

On the other hand, writing 

= h^,hl + Tp 

it is easy to obtain 

2VP- [hlV^h"^) « 2^ + 26^ + 2V''8e^ 

+ y0V^a3p+2V^Tp (All) 

where we have used V^T^ « V^T^ because is first 
order. From Eqs. (AlO, y\ll) wc have 

V'^T^ « -V'^ae^ - 0V^ae^ - V^A^ + ^i? . 

For a sufficiently well behaved i?, which is the case here, 
one can write VPRp = i?, where Rp is actually the gra- 
dient of a scalar. Then 

f ~ ^sip ^ Qsip — Ap + 2^/^ • 
Substituting this back into Eq. (AT), we get 
SuppVPT ^lt{e + V -Si) Hp, - ]^fu,Rp 

-T (apt, + rup^ + V(^ae^) + i/[p£E,] 



(A12) 



-BupV^T. 



(A13) 
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Appendix B: Perturbation Equations in the 
Newtonian Gauge 

The conformal Newtonian gauge can be obtained by 
setting (7 = 0. Defining 



and manipulating Eqs. (55) to (CO), we obtain 

A = 

kz = -3 i^' + n^) , 

T] = -2$. 



(B2) 



With these, and using Eqs. (55) to (GO) and (64) to (08), 
the perturbed field equations in the Newtonian gauge are 
(Bl) derived as 

I 



^nSpf = -frk^^ - 3% ($' + H*) 



1 . 



+ H (f + 2$') + (2-H' + n^)'^ + -fc^(* - *) 

o 



-12^ [n^<i>" + n (3-H' - n^) + -h^*' + {5n' - 2%^) -h^*] 

i^cg/a^ = f/^ - 12/^:^2) k ($' + HvP) + A^I^e-H^^ 



= /Tfc($' + H*) - u^kn in' - H^) $ 
nuf c? = /Tfc^($ - *) + u^kn {n' - n^) x. 



(B3) 



(B4) 
(B5) 

(B6) 
(B7) 



Obviously when /t — 1 = /tt = /ttt = these equa- 
tions reduce to those of general relativity. 
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